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ABSTRACT 



We introduce a new method of searching for and characterizing extra-solar planets. 
We show that by monitoring the center-of-light motion of microlensing alerts using the 
next generation of high precision astrometric instruments the probability of detecting a 
planet orbiting the lens is high. We show that adding astrometric information to the 
photometric microlensing lightcurve greatly helps in determining the planetary mass and 
semi-major axis. We introduce astrometric maps as a new way for calculating astrometric 
motion and planet detection probabilities. Finite source effects are important for low 
Q\ '. mass planets, but even Earth mass planets can give detectable signals. 

Pur 

6 

> 



1 



1. Introduction 

Prompted by several successes the search for extra- 
solar planets has recently intensified. Although the 
radial velocity technique (e.g. Marcy & Butler 1998) 
has yielded the best candidate planets, other promis- 
ing search techniques have been proposed or are un- 



a more complicated "binary lens" caustic structure 
is present and much more complicated photometric 
lightcurves are possible. For a Jupiter-like planet or- 
biting a Sun-like star, the probability of a detectable 
deviation from the standard lightcurve can be sub- 
stantial flMao fc Paczynski 199lfc |Gould fc Loeb 199$ . 



derway: accurate astrometry of nearby stars ( Unwin 



3t al. 1997 ), direct imaging of planets, ( Ftaclas, et al 
1994 ) occultation of stars by their orbiting planets, 
and fine structure on top of photometric microlcnsing 



Encouraged by these calculations, several world-wide 
networks (|Alcock et al. 1997b|; |Albrow et al. 1998[) 



ligfytcurves (Mao & ir'aczynski l99~T; Gould fc l^oeb 



have begun actively searching for such deviations us- 
ing the many dozens of bulge microlensing alerts gen- 
erated yea r ly by the survey experiments ( Alcock, et 



In this paper we propose a new technique that can 
be used to discover and characterize extra-solar plan- 
ets. This method, using astrometric deviations in the 
motion of the center-of- light of microlensing events, 
has several promising features. First, if complete as- 
trometric coverage of microlensing alert events were 
undertaken, we show that the probability of detecting 
any planetary system that is present can be substan- 
tial, typically greater than the probability of detect- 
ing a planetary system when photometric monitoring 
alone is used. Second, if astrometric monitoring is 
added to photometric microlensing we show that in 
most cases the planet mass and projected orbital ra- 
dius can be determined, in many cases with just a 
few astrometric measurem ents. Thus the degenera- 
cies Gaudi fc Gould (1997 ) discussed for photometric 



microlensing can usually be broken by adding astro- 
metric information. This planet detection method, 
like the microlensing photometry method, works best 
for planets at intermediate distances from their star, 

1. e. in the lensing zone. It is sensitive to planets 
throughout the Galaxy, and works well down to quite 
low mass planets. To be useful for this purpose, how- 
ever, the new interferometric instruments and satel- 
lites will need to be able to respond fairly quickly 
to microlensing alerts, and this need may affect their 
design and impact scheduling considerations. 

2. Planet Searching with Photometric Mi- 
crolensing 

Gravitational microlensing is now establ ished as a 
method of detecting low-lumino s ity objects ( Alcock et 
d. 1993j ; [Aubourg, et al. 1993fc [Udalski et al. 1993|) . 
As a low luminosity object passes near the line-of- 
sight to the monitored source, the source brightens in 
a well known achromatic and time-symmetric way. 
If the lens is a faint star with an orbiting planet, 



al. 1997a: Alcock et al. 1996; Udalski et al. 1994 



EROSII 1998|) 



3. Planet Searching with Astrometric Mi- 
crolensing 

Microlensing of a bulge star by a single lens pro- 
duces two images or the star separated by several hun- 
dred microarcseconds and these images move as the 
lens passes near the observer-source line-of-sight. The 
scale of the image motion relative to the source is set 
by the angular Einstein radius r e = 903 /ias[(m; /M Q ) 
(10kpc)(l/D! - 1/Ds)] 1 / 2 , where D s and D t are dis- 
tances to the source and lens respectively. For exam- 
ple, the Einstein radius for a O.3M0 lens at 4 kpc, 
with the source at 8 kpc is 550 /ias. 



Recent work (Gatcwood 1998; Provdo & Shaklan 
1996| ) has shown that the Keck telescope is capable 



of angular resolution down to the sub-milliarcsecond 
scale. In addition, planned interferometers at the 



Keck flColavita, et al. 1998[) and VLT flMariotti, et 



al. 1998 ) should have accuracies of about 10 /zas, 



while NASA's Space Interferome try Mission (SIM) i s 



planned to have 1 /ias resolution ( Unwin, et al. 1997 ) 
These instruments make resolution of astrometric mi- 
crolensing motion a feasible prospect. Several theo- 
retical studies on this subject have been done recently 



(H0g, et al. 1995; 


Miyamoto & Yoshi 1995|; 


Walker 


1995; 


Bodcn, et al. 1998; 


Paczynski 1998; 


Han & 


Chang 1998; 


Han & Tu-Whan 199^). 


Bodcn, et al. 



(1998), for example, describe in detail the predicted 
astrometric motion of lensed images. They point out 
that the instruments described above will only be able 
to measure the motion of the center of light (CoL), 
which is smaller than the motion of the two images, 
but which is nevertheless detectable and useful. The 
CoL motion measured by an unaccelerated observer 
is an ellipse, whose eccentricity is a simple function 
of itmin, the impact parameter of the lens relative to 
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the source. This is depicted in Figure 1. The ellipse, 
however, can be distorted due to the Earth' s motion 
around the Sun (Boden, et al. 1998). Other effects 



such as blending, etc. (Han & Tu-Whan 1998), can 



also distort the ellipse. Figure 1 also shows examples 
of the effects of Earth motion (parallax) and blending 
on the astrometric ellipse. 

Astrometric information is of interest because it 



in only a few hours. As the source crosses a caustic, a 
pair of images is created or destroyed and the magni- 
fication of the these two images becomes very large as 
they approach each other, causing rapid CoL motion. 

Figure ^ shows an example of caustic crossing for 
q = 10 and similar lens parameters as above. This 
corresponds to an earth mass planet around a 0.3M Q 
star. As has been discussed for photometric mi- 



can 


resolve the degeneracies that arise in the pho- 


crolensing (Gaudi & Gould 1997 




Bennett & Rhic 


ton 


tetric microlensing lightcurve. In most cases, a 


1996; 


Wambsganss 1997 


; Gricst & 


Safizadeh 1998), fi- 



fit to the photometric lightcurve gives only the event 
duration which is a function of the three important 
physical quantities: lens mass, distance, and speed. 
Addition of astrometric CoL information, particularly 
in conjunction with astrometric determination of the 
lens parallax, can break this degeneracy and allow 
determination of these physical quantities. 

We note here that a planet orbiting the lens, if po- 
sitioned properly, can perturb the images and thereby 
distort the astrometric ellipse. Figure [2] shows some 
examples of planetary perturbations. We see from 
Figures 1 and || a major difference between planetary 
perturbations and the other distortions - the planet's 
effect is short in duration, while parallax and blend- 
ing are important over the entire duration of the mi- 
crolensing event. This fact enables observers to dis- 
tinguish planetary effects from the other effects. The 
planet's effects are of relatively short duration since 
a planet's Einstein radius is small compared to the 
Einstein radius of its parent star. Although the time 
duration is small, the magnitude of the planetary per- 
turbation can be rather large, even for Earth-mass 
planets, as Figure || illustrates. In Figures |^ and || 
we do not show the parallax and blending effects. We 
assume throughout that they can be independently 
determined. 

Figure ^ shows some examples of planetary astrom- 
etry and photometry curves for q = 10~ 3 (e.g. a Sat- 
urn mass planet orbiting a 0.3M Q star), with D\ = 4 
kpc and D s — 8 kpc. Throughout we will scale to 
an Einstein diameter crossing time of t — 40 days, 



a typical value for bulge microlensing (Alcock, et al 



19£7a). Since all our calculations are done in units 
of the Einstein radius, one can easily rescale them for 
other durations or for other primary lens masses or 
distances. 

Figures ||a and show typical non-caustic cross- 
ing events, with deviations of about 50 /xas lasting 
several days. Figure shows the large rapid motion 
associated with a caustic crossing, here about 200 fia.s 



nite source effects are crucial when considering small 
planets. We have found the same is true for astro- 
metric microlensing. We have calculated finite source 
effects using two different methods. Below we will 
discuss the method of astrometric maps, which au- 
tomatically includes this effect, and which was used 
to calculate the curves in Figure ||. The finite source 
amplifications and centroid motion depicted in Fig- 
ure |3|b, however, were calculated using an ingenious 
technique devised by Gould fc Gaucherel (1997 ), and 
expanded upon by Dominik (1998 ). Figure [3|a shows 
the rapid, large-scale motion of the CoL for a point- 
like source. Note, however, that the planet's effects 
last a far shorter length of time than in Figure ||c. 
Figure ||b shows a close-up of the effects that finite 
sources can have on the astrometric planetary signal. 
Depicted are the astrometric motion for stars of ra- 
dius 1,3,5,9, & 30 i?Q, traversing the same path in the 
source plane as in Figure |^a. As expected, increasing 
the source size smears out the signal, so that for 30 
Rq the deviation is entirely washed out. Notice how- 
ever that for 3 and 5 Rq, the signal is easily visible, 
meaning that Earth mass planets can in principle be 
detected using this technique. 

4. Probability of Planet Detection Using As- 
trometric Microlensing 

The possibility and probability of detecting plan- 
ets with photometric microlensing has been explored 



by several authors (Mao & Paczynski 1991 


; Goulc 


fe Loeb 1992; Wambsganss 1997 


; Bennett & Rhic 


1996; Peale 1997 


; Sackett 1997; 


Sahu 1997 


; Griesi 


k Safizadeh 1998 


) and was found to be substantial. 



We wish to compare astrometric microlensing with 
these studies. In most of these studies, numerous 
lightcurves were generated for a range of planetary 
masses and projected orbital radii, and some simple 
detection criteria was established. In a realistic exper- 
iment, fitting of the planetary lightcurve would need 
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to be performed and the planet mass determined, be- 
fore a planet detection was established. Parameter 
extraction with astrometry will be discussed below, 
but in the following section we will just use simple 
detection criteria analogous to those used in previ- 
ous studies. We define as a detectable perturbation, 
any astrometric curve that contains a deviation from 
the single-lens ellipse greater than a threshold amount 
D t h for a period of time at least tth ■ Detection thresh- 
olds and minimum durations needed will depend upon 
the interferometer or instrument being used, so we 
explore the range of detection thresholds listed in Ta- 
ble 1. For comparison with photometric planetary 
searches, we use a photometric detection c riter ia sim- 
ila r to one used by Bennett fc Rhic (1996 ) and Gricst 



& Safizadeh (1998 ): a photometric deviation of at 



least 4% for a period of at least 6 hours. 

5. Method of Astrometric Maps 

For a planet of mass m p orbiting a lens of mass 
mi, the position of the images z = Xi + iy% are 
found by solving the lens mapping equation z s = 
z + mi / (zi — z) + m p /{z p — z), where Zi, z p , and z s are 
the lens, planetary, and source positions projected to 
the complex lens plane (Witt 199C). For a given z\, 



z p , and z s this equation is a 5th degree polynomial in 
z and has 3 or 5 physical solutions corresponding to 
3 or 5 images. The magnification of each image 
is the reciprocal of the Jacobian determinant of this 
mapping evaluated at the image position z%, and the 
total magnification is just A — ~^2\Ai\. The center- 
of-light (CoL) can be found from zcoL = X) Z »I^M/A 
and it is this quantity that is measured by an interfer- 
ometer. To calculate the CoL motion, one computes 
ZCoL along a given source trajectory through the lens 
plane. This method was used in producing Figures ^ 

D t h (//as) t t h (hours) (t) (hours) 



1 


100 


370 


3 


30 


160 


10 


20 


50 


30 


6 


15 


100 


3 


13 



Table 1: Detection criteria for Figure 3. Listed are 
the deviation threshold D t h, the minimum time ith 
above threshold needed for an event to be counted as 
a detection, and the average time (t) above threshold 
of events passing the cut for x p = 1.3. 



and [|. 

To generate a large number of CoL curves and to 
get an overview of the astrometric motion for all pos- 
sible trajectories, one can alternatively produce an 
"astrometric map" of the source plane using the ray 
tracing technique (Schneider, Ehlers, & Falco 1992, 
p. 303). This is analogous to the magnification maps 
( |Wambsganss 1997] ; priest k Safizadeh 199S| ) pre- 



viously used in calculating photometric microlensing 
probabilities. Using the above lens equation, one sim- 
ply "shoots" many photons from all pixels of the im- 
age back to the source plane, and records the photon 
weighted x and y deflection. The resulting x-shift and 
y-shift maps therefore show the CoL shifts for a source 
at each position in the source planet. A CoL curve 
is created by simply tracing a trajectory through the 
2 maps. To calculate the probability of detecting a 
planet via astrometry, we consider a complete set of 
trajectories (with A max > 1.6) and find the fraction 
of these that pass the thresholds listed in Table 1. 
We then repeat this procedure for various values of 
x p , the projected planet-lens separation, and q, the 
planet-lens mass ratio. We note that since binning of 
the photons is necessary, the finite source-size effect 
is automatically taken into account. Different source 
sizes are easily handled by convolving the maps with 
a round star-sized, limb-darkened kernel. 

The results for a saturn (q = .001) mass planet 
around a 0.3Mq star, with the system at a distance of 
4 kpc, the source star at 8 kpc are shown in Figure |[ 
The astrometric maps used to calculate this figure are 
themselves interesting, and these will be presented in 
a subsequent paper (Safizadeh & Griest 1998). We 



see that the probabilities are substantial for a wide 
range of planet-lens separations. For the case shown, 
x p — 1 corresponds to 2.2 A.U., so the lensing zone 
(0.6 < x p < 1.5) corresponds to 1.3 - 3.3 A.U. 

From Figure [| we see that probabilities range from 
over 70% over the entire lensing zone for the op- 
timistic SIM threshold of 1 /ias deviation down to 
~ 50% for a 3 /Ltas threshold. For the 10 /j,as accuracy 
expected from the Keck and VLT interferometers, the 
probabilities are 20%-40%, while they are below 20% 
for the 100 /xas accuracy reachable with current tech- 
nology. For comparison, the dashed line at around 
30% probability shows detection probabilities using 
photometric lightcurves, and a 4% deviation crite- 
ria. For a jupiter mass planet the probabilities are 
~ 15% higher, while for a 10 earth-mass planet they 
are smaller by a large factor. If the new generation 
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of astrometric instruments do as well as we have as- 
sumed, they should be excellent planet detectors. 

Of course these probabilities should be averaged 
over planetary orbit orientations, and a complete ex- 
ploration of the planetary mass and semi-major axis 
parameter space is necessary. We will present these 
results elsewhere, as well as a study of the the effect 
of the finite source size on small (e.g. Earth mass) 
planets and the loss of signal when the astrometric 
deviation is found by subtracting a fitted astromet- 
ric curve rather than using a priori knowledge of the 
single lens curve. Convolving the saturn mass maps 
with a giant star kernel results in almost unchanged 
probabilities, so finite source effects are not important 
in the saturn mass case presented here. 

6. Extraction of Planetary System Parame- 
ters 

Microlensing has both disadvantages and advan- 
tages with respect to other methods of planet detec- 
tion. The main disadvantage is that the entire signal 
is a short duration deviation on a lightcurve, so fur- 
ther exploration, or even attaining additional infor- 
mation concerning the detected planet, is probably 
impossible. The main advantage of microlensing is 
that large numbers of planets can be found through- 
out the Galaxy and therefore statistics can be found 
on the frequency of planetary systems and the dis- 
tribution of planetary masses and semi-major axes. 
Most other planetary search methods are restricted 
to a small sample of nearby stars. 

Since no follow-up information can be obtained, in 
microlensing it is clearly important to get as much in- 
formation as possible during the brief planetary devia- 
tion from single-source lensing. Thus, adding the two 
additional astrometric curves (Ax and Ay) should be 
extremely valuable. In addition, a convincing demon- 
stration of a microlensing planet detection will require 
determination of the planet mass. This will be done 
via fitting the photometric and/or astrometric curves 
to extract the planetary system parameters. 

To test this hypothesis, we have constructed a sim- 
ple algorithm for fitting data sets. We generated nu- 
merous simulated noisy data sets, and attempted to 
fit them using our automated routine. There are 



8 physical parameters, 5 of which describe the pri- 



parameter it m i n , the Einstein diameter crossing time 
t, the trajectory direction, and the time of closest 
approach. The three planetary parameters are the 
planet-lens mass ratio q, the planet-lens projected 
separation x p , and the planet's angular position rela- 
tive to the lens. 

Since the planet's effects on the astrometric mo- 
tion and magnification are usually perturbative and 
short-lasting, we first fit the data sets to a single iso- 
lated lens, to extract the primary lens parameters. 
We then hold those parameters fixed while we fit for 
the planetary parameters. Given this initial param- 
eter set, we then do a fit in the full 8-dimensional 
parameter space. We find that our automated algo- 
rithm can almost always find a good estimate of the 
planetary parameters, for planetary deviations that 
are significant compared to the added noise. In these 
tests we have ignored planet-lens orbital motion and 
Earth parallax. Orbital motion adds several more pa- 
rameters and will usually be small during the short 
planetary deviation. We note, that if in addition, as- 
trometric parallax ( Bodcn, et al. 1998 ) is found by 
monitoring the event long before or after the plane- 
tary deviation, then the physical size of the lens Ein- 
stein radius can be found and the physical planetary 
mass and separation can be determined. 

Figure |5| shows an example data set and overlaid 
best fit. Note that there are only 8 measurements 
taken during the planetary perturbation, and yet full 
parameter extraction was possible. It is important to 
note that the ability to recover physical parameters 
depends greatly on the noise in the signal. The ex- 
ample of Figure || uses q = 1CP 3 , and x p = 1.3, with 
Gaussian noise of a = 5 /Lias added to each astromet- 
ric data point and 5% Gaussian noise added to each 
photometric data point. If we define a signal/noise 
by the size of the planetary deviation divided by the 
astrometric noise, then we find we are able to recover 
the input parameters reliably for a S/N greater than 
4 or 5 with 6 or 7 data points during the planetary 
deviation. For lower S/ N, we find evidence of the de- 
generacies discussed by Gaudi fc Gould (1997 ) start- 
ing to appear, and for even smaller S/N, parameter 
extraction at all becomes problematic. A more com- 
plete exploration of the effect of noise on the abi lity to 
recover parameters will be presented elsewhere ( Dala| 



Griest I99S 



mary lens and trajectory and three of which describe 
the planet. The 5 primary lens parameters are the 
primary lens Einstein radius, the minimum impact 



Besides noise, the sampling rate can also limit pa- 
rameter extraction. Since time on the new interferom- 
eters will be very valuable and hard to obtain, sam- 
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pling at first may be poor, so it is important to know 
what happens when the assumption of complete as- 
trometric coverage is relaxed. One would like to know 
the minimum number of astrometric measurements 
and the minimum accuracy that is needed to extract 
reasonable values of the planetary parameters. This 
work is in progress. 

As an example of why astrometric data is so help- 
ful, consider the "major image/minor image" degen- 
eracy discussed by paudi fc Gould (1997 ) . For a sin- 
gle lens there are two images. The minor image is 
dimmer and inside the critical curve (Einstein radius) , 
while the major image is always brighter and outside. 
It is not easy to tell from photometry alone whether 
the planet is perturbing the major image or the minor 
image. During a planetary perturbation the bulk of 
the astrometric deviation is due to the changing im- 
age brightness rather than the image motion, so there 
should be a correlation between the photometric and 
astrometric deviations. For example, if the major im- 
age is perturbed and brightens, then the CoL motion 
will be forced outside the single-source ellipse, and the 
photometric deviation will be positive (total magnifi- 
cation increases). However, if the minor image is per- 
turbed and brightens, the CoL motion will be forced 
inside the ellipse with a positive photometric devi- 
ation. Opposite correlations exist for perturbations 
that cause an image to dim, so by simply correlat- 
ing the CoL motion with the photometric deviation 
one easily can tell which image is being perturbed by 
the planet. See Figure || for examples. This corre- 
lation, greatly helps in determining the position of 
the planet, and is used in our algorithm. The above 
correlation works when the perturbation is due to an 
image ncaring the "planetary caustics" , but we find 
we can also distinguish the positio ns for close encoun- 
ters of the "central caustic" (see Griest fc Safizadeh 



(19|98|) for discussion of this terminology). 



The most cri ppling degeneracy discussed by Gaudi 



& Gould (1997 ) arose due to finite source effects. We 
have begun exploring finite source effects, and will 
present our detailed results elsewhere (Dalai & Griest 
1998). For the saturn mass planets of Figure [|, these 
effects are small and our fitting algorithm seems to 
work well. 

7. Summary and Discussion 

By adding astrometric measurements to ongoing 
photometric monitoring of microlensing alerts towards 



the Galactic bulge, one can in many cases, deter- 
mine the planetary mass and projected planet-star 
separation. We feel this is an important reason that 
SIM pointings should be made towards microlens- 
ing alerts. SIM is still in its design phase, and it 
is therefore important that the satellite software and 
hardware be capable of responding quickly to "plane- 
tary alerts" from the photometric microlensing planet 
search teams. We note (from Table 1) that the av- 
erage time above 1 fias deviation for a saturn mass 
planet is about two weeks. On average the time spent 
above 10 ^tas deviation is about two days, while a 30 
/xas deviation lasts 10 hours. For smaller mass planets 
the avera ge durations arc smaller and require faster 
response ( Safizadeh fc Griest 1998 ). 

We have also shown that, using our criteria, we can 
expect to detect and characterize a substantial frac- 
tion of planets that happen to orbit lenses, if reason- 
ably complete astrometric monitoring is undertaken. 
The microlensing surveys typically produce more than 
100 events a year, so microlensing has the potential 
to detect and characterize large numbers of planets. 
Using astrometric and photometric microlensing, we 
could for the first time begin to measure the distribu- 
tion of planet masses and orbital distances throughout 
the Galaxy. 

We thank Andreas Quirrenbach for many helpful 
discussions. This work was supported in part by 
the Department of Energy under grant DEFG0390ER 
40546, and by a Cottrell Scholar award from Research 
Corporation. 
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Fig. 1. — Non-planetary astrometric curves. The 
solid line shows a simple single-lens curve with u m j n = 
0.3, and t — 40 days. The curve is plotted over 
one year, with x's marking each week, so only the 5 
or so weeks at largest y have magnification greater 
than 1.34. The dashed line shows the same with 
an example parallax effect included, while the dot- 
ted line shows the effect of blending (blend fraction 
fb = 60%). 
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Fig. 2. — Some examples of planetary astromet- 
ric and photometric curves. All examples assume 
q = 10 -3 , with a primary lens Einstein radius of 550 
/xas, corresponding to a saturn mass planet. Panel 
(a) has x p = 1.3, panel (b) has x p — 0.7, while panel 
(c) shows a caustic crossing event with x p = 1.3. The 
time axis is scaled so that t — 40 days, and squares 
are plotted one per week, so the durations of the de- 
viations are of order one week. 
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Fig. 3. — Astrometric motion for Earth-mass caustic 
crossing. Panel (a) shows the center-of-light motion 
for a point source, crossing a caustic associated with 
an Earth-mass planet at x p = 0.825. The primary 
lens is 0.3M© at Di — 4 kpc, and source at D s = 8 
kpc. Panel (b) shows a close-up view of the plan- 
etary deviation, with finite-size source. The dotted 
line plots the CoL motion for a 1 Rq size source. The 
solid lines depict the CoL motion for more realistic 
sizes typical of Galactic bulge stars, respectively 3, 5, 
9, & 30 Rq. Note the extreme anisotropy of the axes 
on the graph. For t = 40 days the duration of the 
deviation is about 20 hours, with the center of the 
source spending roughly 90 minutes inside the caus- 
tic. 
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Fig. 4. — Probability of planet detection via astrom- 
etry versus x p , the projected planet-lens separation 
in units of the lens Einstein radius. A planet-lens 
mass ratio q = 1CT 3 and duration (Einstein diame- 
ter crossing time) i = 40 days were assumed. For a 
source distance of 8 kpc, lens distance of 4 kpc, this 
corresponds to a saturn mass planet around a 0.3M Q 
star, or a jupiter mass planet around a 1M star. The 
solid lines are labeled by the thresholds used (detailed 
in Table 1). The dashed line shows the probability for 
photometric detection with a threshold of 4% devia- 
tion over 6 hours. The 1 /ias and 3 /j,as thresholds 
are only relevant for the SIM (Unwin, et al. 1997), 
with 1 yuas being the SIM design goal. The 10 //as, 
and 30 ^tas thresholds would be easy with SIM and 
perhaps possible with the Keck or VLT interferome- 
ters (Colavita, et al. 1998; Mariotti, et al. 1998). The 
100 /ias threshold may be possible with existing in- 
strumentation on the Keck (Gatewood 1998; Provdo 
& Shaklan 1996). 
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Fig. 5. — Recovery of planetary parameters using a 
joint photometric and astrometric fit. The dots show 
simulated data for q = 10~ 3 , x p = 1.4, u m j n = .27, 
with 5 /xas noise added to the astrometric data and 
5% noise added to the photometric data. The solid 
line is the best fit line found by our automated algo- 
rithm, with recovered parameters basically equal to 
those above. Scaling the time axis so t = 40 days, the 
time between each data point is about 9 hours. 



9 



